Numerical atomic orbitals for linear-scaling calculations 
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The performance of basis sets made of numerical atomic orbitals is explored in density-functional 
calculations of solids and molecules. With the aim of optimizing basis quality while maintaining 
strict localization of the orbitals, as needed for linear-scaling calculations, several schemes have 
been tried. The best performance is obtained for the basis sets generated according to a new scheme 
presented here, a flexibilization of previous proposals. The basis sets are tested versus converged 
plane-wave calculations on a significant variety of systems, including covalent, ionic and metallic. 
Satisfactory convergence (deviations significantly smaller than the accuracy of the underlying theory) 
is obtained for reasonably small basis sizes, with a clear improvement over previous schemes. The 
transferability of the obtained basis sets is tested in several cases and it is found to be satisfactory 
as well. 

PACS Numbers: 71.15.Ap, 71. 15. Mb 
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I. INTRODUCTION 

In order to make intelligent use of the increasing power 
of computers for the first-principles simulation of ever 
larger and more complex systems, it is important to de- 
velop and tune linear-scaling methods, where the com- 
putational load scales only linearly with the number of 
atoms in the simulation cell. The present status of these 
methods-and their applications can be found in several 
reviews.Era Essential for linear scaling is locality, and ba- 
sis sets made of localized wavefunctions represent a very 
sensible basis choice. It is not only the scaling what 
matters, however, the prefactor being also important 
for practical calculations. The prefactor depends signifi- 
cantly on two aspects of the basis: (i) the number of basis 
functions per atom, and (ii) the size of the localization 
regions of these functions. 

Atomic orbitals offer efficient basis sets since, even 
though their localization ranges are larger than those 
of some other methods) 3 ] the number of basis functions 
needed is usually quite small. The price to pay for this 
efficiency is the lack of systematics for convergence . Un- 
like with plane-waveo or real-space-gridfl related meth- 
ods, there is no unique way of increasing the size of the 
basis, and the rate of convergence depends on the way 
the basis is enlarged. This fact poses no fundamental 
difficulties, it just means that some effort is needed in 
the preparation of unbiased basis sets, in analogy to the 
extra work required to prepare pseudopotentials to de- 
scribe the effect of core electrons. 

Maximum efficiency is achieved by choosing atomic 
orbitals that allow convergence with small localization 
ranges and few orbitals. It is a challenge again compara- 
ble to the one faced by the pseudopotential-community, 
where transferability and softness are sought .□ For atomic 
wavefunctions the optimization freedom is in the radial 
shape. Gaussian-type orbitals have been proposed for 
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^nnecting with the tradition of quan- 
These bases are, however, quite rigid 
for the mentioned optimization, imposing either many 
gaussians or large localization ranges. 

Numerical atomic orbitals (NAOs) are more flexible 
in this respect. Different ideas have been proposed in 
the literature, originally within tight-binding contexts 
concentrating on minimal (single-^) bases. They are 
obtained by finding the eigenfunctions of the isolated 
atoms confpied-|Within spherical potential wells of differ- 
ent shapes JiSii-a ot-directly modifying the eigenfunctions 
of the free atoms.Ej These schemes give strictly localized 
orbitals, i.e., orbitals that are strictly zero beyond given 
cutoff radii r c . A first extension towards more complete 
basis sets was umposed using the excited states of the 
confined atoms ,E3 but the quite delocalized character of 
many excited states made this approach inefficierrL.unless 
very stringent confinement potentials were used.EZI 

For multiple-C, a better scheme was proposed. -based 
on the split- valence idea of Quantun3_Chemistry,Ejlij but 
adapted to strictly localized NAOs.Ea In the same work, 
a systematic way was proposed to generate polarization 
orbitals suited for these basis sets. The scheme of Rcf. [l8| 
has proven to be quite efficient, systematic, and reason- 
able for a large variety of systems (for short reviews see 
Refs. 

In this work we go beyond previous methodologies be- 
cause of two main reasons: (i) It is always desirable to 
obtain the highest possible precision given the compu- 
tational resources available, and (ii) it is important to 
know and show what is the precision attainable by NAO 
basis sets of reasonable sizes. None of these points was 
systematically addressed in this context before. 

We explore these issues by variationally optimizing ba- 
sis sets for a variety of condensed systems. The param- 
eters defining the orbitals are allowed to vary freely to 
minimize the total energy of these systems. This energy 
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is then compared with that of converged plane-wave cal- 
culations for exactly the same systems, including same 
density functional and pseudopotentials. The optimal 
basis sets are then tested monitoring structural and elas- 
tic properties of the systems. 

The transferability of the basis sets optimized for par- 
ticular systems is then checked by transfering them to 
other systems and testing the same energetical, struc- 
tural and elastic parameters. Finally, the effect of lo- 
calizing the orbitals tighter than what they variationally 
choose is explored on a demanding system. 

II. METHOD 

The calculations presented, bplpw were all done us- 
ing deputy-functional theoryEH (DFT) in its local- 
densityea approximation (LDA). Core electrons were re- 
placed by norm-conseoang pseudopotentialso in their 
fully separable form.ES The— . non-local partial-core 
exchange-correlation correction^ was included for Cu to 
improve the description of the core-valence interactions. 

Periodic boundary conditions were used for all sys- 
tems. Molecules were treated in a supercell scheme al- 
lowing enough empty space between molecules to make 
intermolecular interactions negligible. For solid systems, 
integrations over the Brillouin zone were replaced by con- 
verged sums over selected fc-sets.E3 

Thus far the approximations are exactly the same for 
the two different sets of calculations performed in this 
work: based on NAOs and on plane- waves (PWs). The 
calculations using NAOs were-»pformed with the Siesta 
method, described elsewh.ere.t3E3 Besides the basis itself, 
the only additional approximation with respect to PWs is 
the replacement of some integrals in real space by sums 
in a finite 3D real-space grid, controlled.-by one single 
parameter, the energy cutoff for the gridBj This cutoff, 
which refers to the fineness of the grid, was converged for 
all systems studied here (200 Ry for all except for Si and 
H2, for which 80 Ry and 100 Ry were used, respectively). 
Similarly, the, PW calculations were done for converged 
PW cutoffs.Ell 

Cohesive curves for the solids were obtained by fitting 
calculated energy values for different unit-cell volumes to 
cubic, quartic and Murnaghan-likeo curves, a procedure 
giving values to the lattice parameter, the bulk modulus 
and the cohesive energy of each system. The bulk mod- 
uli given by the Murnaghan and quartic fits deviate from 
each other by around 3%, the Murnaghan values being 
the lowest and the ones shown in the tables. The devia- 
tions between Murnaghan and cubic fits are of the order 
of 7%. The other cohesive parameters do not change 
appreciably with the fits. 

The atomic-energy reference for the cohesive energy 
was taken from the atomic calculations within the same 
DFT and pseudopotentials, always converged in basis 
set. They are hence the same reference for NAOs and 



for PWs, the difference in cohesive energies between the 
two accounting for the difference in the total energy of 
the solid. The isolated-atom calculations included spin 
polarization. 

III. BASIS OF NUMERICAL ATOMIC ORBITALS 

The starting point of the atomic orbitals that conform 
the basis sets used here is the solution of Kohn-Sham's 
Hamiltonian for the isolated pseudo-atoms, solved in a 
radial grid, with the same approximations as for the 
solid or molecule (the same exchange-correlation func- 
tional and pseudopotential). A strict localization of the 
basis functions is ensured either by imposing a boundary 
condition, by adding a confining (divergent) potential, or 
by multiplying the free-atom orbital by a cutting func- 
tion. We describe in the following three main features 
of a basis set of atomic orbitals: size, range, and radial 
shape. 

A. Size: number of orbitals per atom 

Following the nomenclature of Quantum Chemistry, 
we establish a hierarchy of basis sets, from single-^ to 
multiple-^ with polarization and diffuse orbitals, cover- 
ing from quick calculations of low quality to high preci- 
sion, as high as the finest obtained in Quantum Chem- 
istry. A single-^ (also called minimal) basis set (SZ in 
the following) has one single radial function per angular 
momentum channel, and only for those angular momenta 
with substantial electronic population in the valence of 
the free atom. It offers quick calculations and some in- 
sight on qualitative trends in the chemical bonding and 
other properties. It remains too rigid, however, for more 
quantitative calculations requiring both radial and angu- 
lar flexibilization. 

Starting by the radial flexibilization of SZ, a better 
basis is obtained by adding a second function per chan- 
nel: double-^ (DZ). Several schemes have been proposed 
to generate this seconds-function. In Quantum Chem- 
istry, the split valencetQtB scheme is widely used: start- 
ing from the expansion in Gaussians of one atomic or- 
bital, the most contracted gaussians are used to define 
the first orbital of the double-^ and the most extended 
ones for the second. Another proposal defines the sec- 
ond £ as the, derivative of the first one with respect to 
occupation.E3_For strictly localized functions there was a 
first proposaO of using the excited states of the confined 
atoms, but it would work only for tight confinement. An 
extension of the split valence idea of Quantum Chemistry 
to strictly localized NAOs was proposed in Ref. |l^ and 
has been used quite succesfully in a variety of systems. 
We follow this scheme in this work, which generalizes to 
multiple-^ trivially. 
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Angular flexibility is obtained by adding shells of 
higher angular momentum. Ways to generate these so- 
called polarization orbitals have, been described iu_the 
literature, both for GaussiansEjO and for NAOs.Ea In 
this work, however, they will be obtained variationally, 
as the rest, within the flexibilities described below. 

B. Range: cutoff radii of orbitals 

Strictly localized orbitals (zero beyond a cutoff radius) 
are used in order to obtain sparse Hamiltonian and over- 
lap matrices for linear scaling. The traditional alterna- 
tive to this is based on neglecting interactions when they 
fall below a tolerance or when the atoms are beyond some 
scope of neighbours. Its disadvantage is that it implies 
a deviation from the original Hilbert space, the total en- 
ergy is no longer a variational magnitude, and produces 
numerical instabilities. 

For the bases made of strictly localized orbitals, the 
problem is finding a balanced and systematic way of 
defining all the different radii, since both the precision 
and the computational efficiency in the calculations de- 
pend on them. A scheme was proposedES in which all 
radii were defined by one single parameter, the energy 
shift, i.e., the energy raise suffered by the orbital when 
confined. In this work, however, we step back from that 
systematic approach and allow the cutoff radii to vary 
freely in the optimization procedure (up to a maximum 
value of 8 a.u.). 

C. Shape 

Within the pseudopotential framework it is important 
to keep the consistency between the pseudopotential and 
the form of the pseudoatomic orbitals in the core region. 
This is done by using as basis orbitals the solutions of the 
same pseudopotential in the free atom. The shape of the 
orbitals at larger radii depends on the cutoff radius (see 
above) and op,the way the localization is enforced. The 
first proposalO used an infinite square-well potential. It 
has been widely and succesfully used for minimal bases 
within the ab-initio tight-binding scheme of Sankey and 
collaboratorsE£l using the Fireball program, but also 
for more flexible bases using the methodology of Siesta. 

This scheme has the disadavantage, however, of gener- 
ating orbitals with a discontinuous derivative at r c . This 
discontinuity is more pronounced for smaller r c 's and 
tends to disappear for long enough values of this cutoff. 
It does remain, however, appreciable for sensible values 
of r c for those orbitals that would be very wide in the 
free atom. It is surprising how small an effect such kink 
produces in the total energy of condensed systems (see 
below). It is, on the other hand, a problem for forces and 
stresses, especially if they are calculated using a (coarse) 
finite three-dimensional grid. 



Another problem of this scheme is related to its defin- 
ing the basis considering the free atoms. Free atoms can 
present extremely extended orbitals, their extension be- 
ing, besides problematic, of no practical use for the calcu- 
lation in condensed systems: the electrons far away from 
the atom can be described by the basis functions of other 
atoms. 

Both problems can be addressed simultaneously by 
adding a soft confinement potential to the atomic Hamil- 
tonian used to generate the basis orbitals: it smoothens 
the kink and contracts the orbital as variationally suited. 
Two soft confinement potentials have been proposed in 
the literature, both of the form V(r) = V r n , one for 
n=2 (Ref. |l|) and the other for n=6 (Ref. 0). They 
present their own inconveniences, however. Firstly, there 
is no radius at which the orbitals become strictly zero, 
they have to be neglected at some point. Secondly, these 
confinement potentials affect the core region spoiling its 
adaptation to the pseudopotential. 

This last problem affects a more traditional scheme as 
well, namely, the one based on the radial scaling of the 
orbitals by suitable scaie, factors. In addition to very ba- 
sic bonding arguments, Eil it is soundly based on restoring 
virial's theorem for finite bases, in the case of coulombic 
potentials (all-electron calculations )£3 The pseudopoten- 
tials limit its applicability, allowing only for extremely 
small deviations from unity (~ 1%) in the scale factors 
obtained variationally (with the exception of hydrogen 
that can contract up to 25%) H ._. 

An alternative scheme has also been proposed:E3 In- 
stead of modifying the potential, it directly modifies the 
orbitals of the free atom. Following ideas of previous 
mixed-basis schemesE3 the atomic orbital is multiplied 
by 1 — exp[— a(r — r c ) 2 ] for r < r c and zero otherwise. 
This scheme does provide strict localization beyond r c , 
but introduces a different problem: for large a and small 
r c a bump appears in the orbital close to r c , which be- 
comes a discontinuity in the wave-function in the limit of 
zero a. 

In this work we propose a new soft confinement poten- 
tial avoiding the mentioned deficiencies. It is flat (zero) 
in the core region, starts off at some internal radius rj 
with all derivatives continuous and diverges at r c ensur- 
ing the strict localization there. It is 



V{r) = V o e —^- (1) 
r c - r 

In the following the different schemes are compared, their 
defining parameters being allowed to change variation- 
ally. 

Finally, the shape of an orbital is also changed by the 
ionic character of the atom. Orbitals in cations tend to 
shrink, and they swell in anions. Introducing a 5Q in 
the basis-generating free-atom calculations gives orbitals 
better adapted to ionic situations in the condensed sys- 
tems. 
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TABLE I. Comparison of different confinement schemes on 
the cohesive properties of MgO, for SZ and DZP basis sets. 
The generalization of the different schemes to DZP is done as 
explained in the text. Unconfined refers to using the uncon- 
fined pseudoatomic orbitals as basis, a, B, and E c stand for 
lattice parameter, bulk modulus, and cohesive energy, respec- 
tively. The PW calculations were performed with identical 
approximations as the NAO ones except for the basis. Aug- 
mented plane-wave (LAPW) results were taken from Ref. UA, 
and the experimental values from Ref. 43l 
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IV. OPTIMIZATION PROCEDURE 

Given a system and a basis size, the range and shape 
of the orbitals are defined as described above, depend- 
ing on parameters. These are defined variationally: the 
energy is minimized with respect to all of them. As a ro- 
bust and simple minimization method not requiring the 
evaluation of derivatives, we have chosen the downhill 
simplex method.Ell 

The parameters varied in an optimization are the fol- 
lowing. Per atomic species there is a global 5Q. Confine- 
ment is imposed differently for each angular momentum 
shell, with its corresponding parameters that depend on 
the scheme used. Hard confinement implies one parame- 
ter per shell (r c ), and our soft confinement implies three 



(r c , r.i, and V Q ). One para: 




(V ) is needed only in 
the r n -confinement schemesJIS'Eil and two parameters in 
the scheme of Kenny et al£3 (r c and the width of the 
cutting function). Finally, for each £ beyond the first, 
there is a matching radius as defined in Ref. The 
values obtained for the parameters in the optimizations 
described below can be obtained from the authors. 



V. RESULTS 

A. Comparison of different confinement schemes 

Table | shows the performance for MgO of the differ- 
ent schemes described above for constructing localized 
atomic orbitals. The basis sets of both magnesium and 
oxygen were variationally optimized for all the schemes. 
Mg was chosen because the 3s orbital is very extended in 
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Free atom 

This work 
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FIG. 1. Shape of the 3s orbital of Mg in MgO for the differ- 
ent confinement schemes (a) and corresponding potentials (b). 



the atom and both the kink and the confinement effects 
due to other orbitals are very pronounced. Results are 
shown for a SZ (single s and p channels for both species) 
and a DZP basis (double s and p channels plus a single d 
channel). Fig. [j] shows the shape of the optimal 3s orbital 
for the different schemes, and the shape of the confining 
potentials. 

The following conclusions can be drawn from the re- 
sults: (i) Within the variational freedom offered here, 
the 3s orbital of Mg wants to be confined to a radius 
of around 6.5 Bohr, irrespective of scheme, which is ex- 
tremely short for the free atom. This confinement pro- 
duces a pronounced kink in the hard scheme, (ii) The 
total energy is relatively insensitive to the scheme used 
to generate the basis orbitals, as long as there is effective 
confinement. (Hi) The basis made of unconfined atomic 
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TABLE II. Basis comparisons for bulk Si. a, B, and E c 
stand for lattice parameter (in A), bulk modulus (in GPa), 
and cohesive energy (in eV), respectively. SZ, DZ, and TZ 
stand for single-^, double-^, and triple-^". P stands for po- 
larized, DP for doubly polarized. LAPW results were taken 
from Ref. kg, and the experimental values from Ref. by. 
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orbitals is substantially worse than any of the others, (iv) 
The pronounced kink obtained in Sankey's hard confine- 
ment scheme is not substantially affecting the total en- 
ergy as compared with the other schemes. It does per- 
turb, however, by introducing inconvenient noise in the 
energy variation with volume and other external param- 
eters, and especially in the derivatives of the energy, (v) 
The scheme proposed in this work is variationally slightly 
better than the other ones, but not significantly. Its main 
advantage is the avoidance of known problems. In the re- 
mainder of the paper, the confinement proposed in this 
work will be used unless otherwise specified. 



B. Basis convergence 

Table fi| shows how NAO bases converge for bulk sil- 
icon. This is done by comparing different basis sizes, 
each of them optimized. The results are compared to 
converged (50 Ry) PW results (converged basis limit) 
keeping the rest of the calculation identical. They are 
also compared to all-electron LDA resultsCj to compare 
basis errors with the ones produced by the pseudopoten- 
tials. Experiment gives a reference to the accuracy of the 
underlying LDA. Fig. ||(a) shows the cohesion curves for 
this system. 

The comparisons above are made with respect to the 
converged-basis limit, for which we used PWs up to very 
high cutoffs. It is important to distinguish this limit from 
the PW calculations at lower cutoffs, as used in many 
computations. To illustrate this point, Fig. ||(b) com- 
pares the energy convergence for PWs and for NAOs. 
Even though the convergence of NAO results is a priori 
not systematic with the way the basis is enlarged, the 
sequence of bases presented in the Figure shows a nice 
convergence of total energy with respect to basis size (the 
number of basis functions per atom are shown in paren- 
thesis in the Figure): the convergence rate is similar to 
the one of PWs (DZP has three times more orbitals than 
SZ, and a similar factor is found for their equivalents in 
PWs). For the particular case of Si, Fig. |^ shows that 
the polarization orbitals (3d shell) are very important for 
convergence, more than the doubling of the basis. This 
fact is observed from the stabilization of SZP with respect 
to SZ, which is much larger than for DZ. 
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(25) (71) (130) (201) (280) — PW basis size 
FIG. 2. Convergence of NAO basis sets for bulk Si. 

(a) Cohesive curve for different basis sets. The lowest curve 
shows the PW results, filled symbols the NAO bases of 
this work (opt), and open symbols the NAO bases following 
Ref. |l8|. Basis labels are like in Table [o| F standing for an 
extra /-shell. 

(b) Comparison of NAO convergence with PW convergence. 
In parenthesis the number of basis functions per atom. 
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TABLE III. Equivalent PW cutoff (E cut ) to optimal DZP 
bases for different systems. Comparison of number of basis 
functions per atom for both bases. For the molecules, a cubic 



unit cell of 10 A of side was used. 

System # funct. DZP # funct. PW -E cut (Ry) 

H2 5 11296 34 

2 13 45442 86 

Si 13 144 16 

Diamond 13 284 59 

a-quartz 13 923 76 



Fig. g shows that an atomic basis requires ten times less 
functions than its (energetically) equivalent PW basis, 
being Si the easiest system for PWs. For other systems 



TABLE IV. Basis comparisons for different solids, a, B, 
and E c stand for lattice parameter (in A), bulk modulus (in 
GPa), and cohesive energy (in eV), respectively. 



the ratio is much larger, as shown in Table III. 

It is important to stress that deviations smaller than 
the ones due to the pseudopotential or the DFT used 
are obtained with a relativelymodest basis size as DZP. 
This fact is clear in Table 10 for Si, and in Table IV 
for other systems. Table [V summarizes the cohesion 



results for a variety of solids of different chemical kind. 
They are obtained with optimal DZP basis sets. It can 
be observed that DZP offers results in good agreement 
with converged-basis numbers, showing the convergence 
of properties other than the total energy. The deviations 
are similar or smaller than those introduced by LDA or 
by the pseudopotentialo 



VI. TRANSFERABILITY 

To what extent do optimal bases keep their perfor- 
mance when transfered to different systems than the ones 
they were optimized for? This is an important question, 
since if the performance does not suffer significantly, one 
can hope to tabulate basis sets per species, to be used for 
whatever system. If the transferability is not satisfactory, 
a new basis set should then be obtained variationally for 
each system to be studied. Of course the transferability 
increases with basis size, since the basis has more flexibil- 
ity to adapt to different environments. A thorough study 
of transferability, therefore, implies many calculations in 
different systems, for different bases and different basis 
sizes. In this work we limit ourselves to try it on DZP 
bases for a few representative systems. 

Satisfactory transferability has been obtained when 
checking in MgO the basis set optimized for Mg bulk and 
O in a water molecule. Similarly the basis for O has been 
tested in H2O and O2, and the basis for C in graphite 
and diamond. Again, the results show deviations due to 
the basis that are smaller than the errors introduced by 
the pseudopotentials and/or the DFT functional. The 
results are shown in Table |v|. 

Table VI shows the results for the structural parame- 
ters of Si02 in its a-quartz structure. The DZP numbers 
have been obtained for a basis that was optimized not 
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a- quartz. 



Tightening the confinement of the basis in 



TABLE V. Transferability of basis sets. "Transf" stands 
for the DZP basis transferee! from other systems, while "Opt" 
refers to the DZP basis optimized for the particular system. 
For MgO the basis was transfered from bulk Mg and an H2O 
molecule, for graphite the basis was transfered from diamond, 
and for H2O it was taken from H2 and O2. a, B, and E c 
stand for lattice parameter, bulk modulus, and cohesive en- 
ergy, respectively. AE stands for the energy difference per 
atom between graphite and a graphene plane. Et is the bind- 
ing energy of the molecule. 



System 


Basis 


Properties 


MgO 




a (A) 


B (GPa) 


E c (eV) 




Transf 


4.13 


157 


11.81 




Opt 


4.10 


167 


11.87 




PW 


4.10 


164 


11.89 




LAPW 


4.26 


147 


10.40 




Exp 


4.21 


152 


10.30 


Graphite 




a (A) 


c(A) 


AE (meV) 




Transf 


2.456 


6.50 


38 




PW a 


2.457 


6.72 


24 




Exp b 


2.456 


6.674 


23 c 


H 2 




d -H (A) 


6>o-h-o (deg) 


E b (eV) 




Transf 


0.975 


105.0 


12.73 




Opt 


0.972 


104.5 


12.94 




PW 


0.967 


105.1 


13.10 




LAPW d 


0.968 


103.9 


11.05 




Exp 6 


0.958 


104.5 


10.08 



a M. C. Schabel, |Ref. 
b Y. Baskin, RefEi 
C L. A. GirifalcoJjtef. 
d P. Serena, RefB p 
C G. Herzberg, Ref.E3 



E3 



TABLE VI. Performance of the basis of Si and O as op- 
timized in c-Si and in a water molecule, respectively, for the 
structural parameters of a-quartz. 







Exp a 


PW b 


PW C 


PW d 


PW C 


DZP 


a (A) 




4.92 


4.84 


4.89 


4.81 


4.88 


4.85 


c(k) 




5.41 


5.41 


5.38 


5.32 


5.40 


5.38 


^Si-O 


(A) 


1.605 


1.611 


1.60 


1.605 




1.611 


d 2 

"Si-O 


(A) 


1.614 


1.617 


1.60 


1.605 




1.611 


a?si-o 


-Si (deg) 


143.7 


140.2 




139.0 




140.0 



a L. Levien, C. T. Prewitt and D. J. Weidner, Ref. || p 
b P. Sautet (unpublished) using ultrasoft pseudopotentials.Ej 
C D. R. Hamann, Ref. |52| . 
d G.-M.Rignanese et aL, Ref. ||. 

C F. Liu et aL, Ref. BOl using ultrasoft pseudopotentials. 



s 


(a.u.) 




To 


(a.u 


■) 


a 


c 


,1 

"Si-O 


,2 

"Si-O 


OSi-O-Si 


S 


P 


d 


S 


P 


d 


(A) 


(A) 


(A) 


(A) 


(deg) 


8.0 


8.0 


8.0 


8.0 


8.0 


8.0 


4.85 


5.38 


1.611 


1.612 


140.0 


6.0 


6.0 


6.0 


8.0 


8.0 


8.0 


4.85 


5.35 


1.607 


1.608 


140.0 


6.0 


6.0 


6.0 


5.0 


5.0 


5.0 


4.74 


5.29 


1.610 


1.610 


134.0 


6.0 


6.0 


6.0 


4.5 


4.5 


4.5 


4.69 


5.26 


1.610 


1.610 


132.0 


6.0 


6.0 


6.0 


5.0 


6.5 


4.0 


4.84 


5.36 


1.607 


1.608 


139.7 


5.6 


6.3 


4.2 


4.0 


5.3 


2.8 


4.81 


5.34 


1.607 


1.610 


138.2 



for ev-quartz itself, but for bulk silicon for the Si basis 
and for the water molecule for the O basis. SiC>2 was 
chosen because its being quite sensitive to many approx- 
imations and in particular to the basis sat. It was hard 
to converge for previous NAO schemes ,c3 givingcj typi- 
cally longer Si-O bonds (with deviations of around 1.5%) 
and smaller unit cells (deviations of around 1.5% and 2% 
for the a and c parameters, respectively). The results of 



Table VI are very satisfactory, showing (i) the good per- 
formance of NAOs, (ii) their transferability in this case, 
and (Hi) the improvement of the basis sets proposed here 
over previous bases. 



VII. LIMITING THE RANGE 

In this work we have concentrated on variationally op- 
timized basis sets, allowing the cutoff radii for the differ- 
ent orbitals to vary freely, as long as the orbitals remained 
strictly localized. This was done in the spirit of exploring 
the capabilities of the NAO basis sets. Some orbitals de- 
manded reasonably short values of r c , others chose long 
ranges. As mentioned earlier, the range of the orbitals is 
important for the efficiency in the calculations. There- 
fore, further work will be very important to explore the 
possibility of enforcing smaller ranges in reasonably bal- 
anced ways and its effect on the precision. A systematic 
study in this direction will be subject of future work, we 
have limited ourselves here to illustrate the nature of the 
problem in the particular example of a-quartz. 

The basis has been optimized as before (Si in bulk Si 
and O in H2O), but imposing now tighter r c 's. The re- 
sults are summarized in Table VII . The constraining of 
Si orbitals to 6.0 a.u. affects the geometry only slightly, 
whilst the contraction of the O orbitals to 5.0 a.u. im- 



plies a substantial contraction of the cell due to the de- 
crease of Si-O-Si angle rather than the shortening of the 
Si-O bond. Note that, from an atomic perspective, the 
confinement of Si to 6 a.u. is tighter than the 5 a.u. 
confinement of O. 

Allowing for different r c 's for the different channels we 
observe that the shrinking of the cell is avoided keeping 
a long p orbital for O, the s and cLremaining comparably 
shorter, indicating the (expectecliS) different "compress- 
ibilities" of the different orbitals. As a candidate of unify- 
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ing criterion, we tested the simple definition of r c for each 
channel as the value of the radius at which the uncon- 
strained orbital would have a defined value (0.01 a.u.~ 3 ). 
Even though the energy raise was appreciable (about 70 
meV per atom) , the geometry retained an acceptable pre- 
cision, the orbitals being quite short, thus allowing quite 
efficient calculations. Further work is, however, needed 
to explore in detail this and other possibilities. 



VIII. CONCLUSIONS 

The variational optimization of NAO basis sets for 
different systems allows us to draw the following con- 
clusions, (i) The performance of NAO basis sets of 
modest size as DZP is very satisfactory for the systems 
tried. For this basis size, the deviation from basis conver- 
gence is smaller than errors due to the pseudopotential 
or DFT used, (ii) The bases obtained here represent a 
substantial improvement over previous NAO basis sets. 
In particular, the optimization in condensed systems of- 
fers better and more efficient bases than purely atomic 
schemes. (Hi) The radial shapes of the orbitals obtained 
as proposed in this work offer better bases than previ- 
ous schemes from a variational point of view, albeit not 
a substantial difference is obtained, (iv) The elimina- 
tion of the discontinuity in the derivative, while retaining 
strict localization and leaving the core region untouched, 
gives bases of better quality from the point of view of 
the energy, its derivatives, and computational efficiency. 
(v) The bases obtained showed enough transferability to 
expect that a basis tabulation would be useful, and that 
the optimization of the basis for each particular system 
will not be necessary. Finally, the selective sensitivity 
to orbital-range tightening has been shown, making clear 
the need of further work systematically to control the 
cutoff radii for improving efficiency without loss of pre- 
cision. 
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